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Abstract 

Experiments on charm hadroproduction have shown a substantial 

^^ , difference in the production of charm and anticharm hadrons. In this 

work we study A+ and A 7 inclusive production in p— N interactions in 

r S^ ■ the framework of two component models. We show that the recombi- 

^h, nation two component model gives a qualitatively and quantitatively 

Q_i! good description of Ajf production while the intrinsic charm model 

seems to be ruled out by recent experimental data from the SELEX 

Collaboration. 
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It is well known from experiments that there is a substantial difference 
in the production of charm and anticharm hadrons in hadron-hadron inter- 
actions. Leading (L) particles, which share one or more valence quarks with 
the initial hadrons, are favored in the incident hadron direction over Non- 
Leading (NL) particles, which share none. This effect, known as Leading 
Particle Effect, has been observed in inclusive production of charm mesons 
and baryons in tt~ — N, p — N, K — N and S~ — N interactions by several 
experiments [|I|, 0, ||] . 

Leading particle effects can be quantified by means of a production asym- 
metry, which is denned as 

" da L + da NL " l ' 
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Asymmetries in charm-anticharm production indicate that charm hadro- 
nization cannot proceed by independent fragmentation alone. They imply 
also that some sort of recombination mechanism, involving valence quarks in 
the initial hadrons, must take place in the production. 

Several models have been proposed to explain charm hadropro duct ion. 
However, no theoretical consensus has been reached yet on what is the dom- 
inant mechanism giving rise to the observed asymmetries. This is partly due 
to the lack of simultaneous measurements of charm-anticharm asymmetries 
and inclusive particle distributions in the same experiment. Actually, as 
models depend on a set of parameters, and these parameters can be adjusted 
to describe either the charm-anticharm asymmetries or the inclusive particle 
distributions, a meaningful comparison of models to experimental data must 
be done on both, asymmetries and inclusive particle distributions. 

Among the proposed models are the String Fragmentation model (SF) M, 
implemented in the Lund Pythia-Jetset package ||; the recombination of 
charm quarks produced perturbatively with the remnants of the initial ha- 
drons |J; the Intrinsic Charm model (IC) J7| and the recombination two 
component model (R2C) |§. 

Aiming to extract information on the charm hadron production mecha- 
nisms, in this work we shall compare predictions of both, the IC and R2C 
two component models to recent experimental data on p + N — > Af + X by 
the SELEX Collaboration §. 

In two component models the total cross section for charm hadron pro- 
duction receives contributions from two different processes, namely perturba- 
tive production of a cc pair in QCD followed by independent fragmentation, 
and contributions coming from another, non-perturbative, mechanism. So 
far, two possibilities have been considered for this second contribution: IC 
coalescence |7|] and the recombination of charm quarks, already present in 
the sea of the initial hadrons, with valence and sea quarks from the initial 
hadrons |J . Then, in the p + N — > Af + X reaction, 

da da Fra 9- da IC( - Rec ^ 
= + ■ ( 2 ) 

(JjJbp CLJb p (XJb ]7 

The first term in the RHS of Eq. (|2]) describes the production of charm 
hadrons by independent fragmentation, while the second one is the contri- 
bution to A^ inclusive production coming from the IC or the recombination 
mechanisms. 



Charm hadron production by independent fragmentation is given by 
da Fraa - 1 



dx F 2 



^jH ab (x a ,x b ,Q 2 ) 



where 



x LE±d*idzdi%dy , (3) 
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DkA*) = n : , — /n ^ (4) 
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is the Peterson fragmentation function || with N\ c a normalization constant. 
H ab (x a , Xb, Q 2 ) contains information on the initial hadron structures and the 
dynamics of the perturbative QCD (pQCD) charm production. At Leading 
Order (LO) it is given by 

H ab (x a , x b ,Q 2 ) = £ a , b q a (x a ,Q 2 )q b (x b ,Q 2 ) 

da , 



+ q a (x a ,Q )qb(xb,Q 



dt m 
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+ 9a(x a ,Q 2 )g b (x b , Q^-jlgg . (5) 

dt 

In Eqs. (0) (|) and @, Xf, i = a,b; is the momentum fraction of the initial 
hadron carried by the parton i, z and p 2 ^ are the momentum fraction of 
the initial hadron carried by the charm quark and its transverse momentum 
squared respectively, and y is the rapidity of the charm antiquark. The sum 
in Eq. (|5|) runs over light and strange quarks. 

Up to LO in the pQCD processes (see Fig. 0), no contribution to the 
A+ — A^ asymmetry arises from the charm quark production. At Next to 
Leading Order (NLO), a small c — c asymmetry translates into a tiny A+ — A~ 
asymmetry |TTJ. However, this effect is very small and has the opposite sign 



to the experimentally observed asymmetry in p + N interactions. Then, the 
second term in the RHS of Eq. (Q) must give the dominant contribution to 
the A+ — A~ asymmetry. 

In the IC model for the p+N -^ Af+X reaction in the x F > region, the 
second term in the RHS of Eq. (0) comes from fluctuations of protons in the 
beam to \uudcc) Fock states [j/j. These Fock states break up in the collision 
contributing to A+ production through the coalescence of the intrinsic charm 
quark with u and d quarks. To obtain a A^ purely from this process, a 
fluctuation of the proton to a uuduuddcc) Fock state is required. As the 






Figure 1: Parton fusion processes contributing to cc perturbative production 
at LO. 



probability of the later is smaller than for the former, A+ production is 
favored in the proton direction. The A+ differential cross section for the 
intrinsic charm process is |7| 



where 



da IC 



dx. 



dP IC 






X yX F X u Xd X c ) 
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dp IC 



X 



N 5 at (M c 2 5 ) 
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(6) 



(7) 



is the probability of the \uudcc) fluctuation of the proton, and r is a param- 
eter which must be fixed from experimental data. Neglecting contributions 



coming from the uuduuddcc) Fock states of the proton, which are very small, 
the A~ differential cross section is given only by the first term in the RHS of 
Eq. (i. 

In the R2C model, the second term in the RHS of of Eq. (fj) has the 
form H |TT 



CMJ j^ ec j CLX U CLX d CLX C . , 
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x R 3 (x u ,x d ,x c ,x F ) , (8) 

where F% and -R3 are the multiquark distribution and the recombination 
functions respectively. As for the intrinsic charm model, r Rec is a parameter 
which must be fixed from experimental data. 

For A+ production in p + N interactions in the x F > region, F3 is given 
by 

F 3 (x u , x d , x c ) ~ x u u(x u )x d d(x d )x c c(xc) 

x p(x u ,x d ,x c ) , (9) 

while for A~ production 

-^3 \ x ui % d i Xc) ~ XuU(Xu)x d d{X d )x 5 c(x 5 ) 

X p(xu,x d ,x E ) ■ (10) 

In Eqs. (|j) and ([10|), xq(x) is the q-flavored quark distribution in the proton. 
Notice that the multiquark distribution of Eq. @ receives contributions from 
u and d valence and sea quarks, while Eq. (|H]) is constructed only by using sea 
quark distributions. Thus, A~ production is due solely to the recombination 
of antiquarks popped up from the vacuum in the interaction. As quarks and 
antiquarks are created in pairs from the vacuum, A+'s can also be formed by 
the recombination of u and d sea quarks, in addition to the u and d valence 
quarks, with charm quarks, thus giving rise to a A+ — A~ asymmetry. 

The momentum correlation function p, which we assume to be the same 
for both A+ and A~ production, is usualy taken as |TlJ 



p(x u ,x d ,x c ) = (1 - x u -x d -x c y , (11) 

with the exponent 7 fixed appealing to some consistency condition like fTT| 
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for the valence quarks in the proton. 

In Eqs. ([|) and (|To|) we are assuming the existence of charm quarks inside 
the proton. Indeed, assuming that the global scale of the whole process is of 
the order of Q 2 ~ 4m);, there should be a substantial contribution of charm 
quarks to the proton structure ||12|| . However, charm quarks in the proton 
may have a twofold origin, namely, a non-perturbative component which must 
exist over a time scale independent of Q 2 , and a perturbative component due 
to the QCD evolution. The non-perturbative contribution is expected to be 
small, of the order of 1 % or less [|H|, and at Q 2 ~ Am 2 the perturbative 
component must be dominant. On the other hand, by assuming the existence 
of charm inside the proton we are consistently including the flavor exitation 
diagrams which are not considered in the LO calculation of Eq. (fj)- ([5]) |14[ 
(See Fig. ||]). Note that flavor exitation diagrams are usually included in 
the pQCD calculation at NLO, but only for the perturbatively generated 
charm quarks. The non-perturbative charm sea must be taken into account 
through the recombination process. Moreover, it is difficult to account for 
the recombination of the spectator c-quark (see Fig. |2|) when flavor exitation 
diagrams are included into a NLO pQCD calculation of charm production. 

The recombination function R3 is given by | 1~5| 
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with rii — 1 and n-i = 5 JT6| and a is a normalization constant. The same 
recombination function is used for A~ inclusive production. 

In order to compare model predictions to experimental data, we have 
used 



dN A < 
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da 
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and similarly for the A~ differential cross section. In the equation above, 
N is a global normalization constant which has been adequately fixed from 
experimental data. The r IC,Rec parameter was allowed to be different for the 
IC and R2C models, however, it was fixed to the same value for particle and 
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Figure 2: Upper: A+ production by recombination of a c quark in the sea of 
the proton with u and d quarks. Lower: flavor exitation diagrams contribu- 
ting to A+ and A~ production by recombination. 

antiparticle production within each model. The A+ — A~ asymmetry as a 
function of x F was calculated according to Eq. (fj) using the parametrization 
of Eq. fll4|) for the particle and antiparticle x F inclusive distributions. 

To control the size of each individual contribution to the total differential 
cross section, the calculated distributions for the fragmentation, IC and re- 
combination processes for A+ production were each normalized to unity. The 
recombination differential cross section for A~ production was normalized by 

multiplying by a -1 , where a = J -^- c dx before normalization. Thus, the 
contribution of the light sea to the A+ and A~ production is the same. In 
this way, the coefficients r IC and r Rec give the relative size of the IC and 
Recombination component respectively, in comparison to the independent 
fragmentation process in the total cross section. 

Prediction by the IC and R2C models are shown in Fig. |3] for both the 
A+ — A^ asymmetry and the inclusive particle distribution as a function of 



x F . Model predictions are also compared to the p + N 
the SELEX Collaboration HI. 



A^ + X data from 



Curves were obtained using the GRV-92 |T2[ parton distributions in nu- 
cleons in Eqs. (|5|), (|9|) and fllPl). The exponent 7 in the correlation function 
of Eq. (|TT[) was fixed to 7 = —0.1 || and the overall Q 2 scale was chosen as 
Q 2 = 4m 2 with m c = 1.5 GeV. In the Peterson fragmentation function we 
used e = 0.06. In order to fix the parameter r IC,Rec in both the IC and R2C 
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Figure 3: Left: A+ — A 7 asymmetry in pN interactions. Solid line is the 
prediction of the R2C model with r Rec = 1.5 and dashed line is the prediction 
of the IC model with r = 0.4. Right: A c xf distribution in pN interactions. 
Solid line is the prediction of the R2C model and dashed line is the prediction 
of the IC. Dot-dashed line shows the A~ distribution as predicted by the R2C 
model. Experimental data are from Ref. [3]. 

models, we adjusted the curves to describe first the asymmetry. Once this 
parameter was fixed, we compared model predictions to experimental data 
on the x F particle distribution. 

As can be seen in the figure, the IC model cannot describe simultaneously 
the A+ — A~ asymmetry and the A c x F distribution. On the other hand, the 
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R2C model gives a qualitative and quantitatively good description of both, 
the asymmetry and the x F particle distribution. 

Furthermore, in order to fit the A c x F distribution, the recombination 
part of Eq. (|2|) must be of the order of 1.5 times bigger than the fragmenta- 
tion contribution. This implies that the main mechanism in A c production is 
recombination, even in the low x F region (see Fig. ^). It is interesting to note 
also that recombination seems to be bigger than independent fragmentation 
also for A~ production. The IC contribution to the proton structure, which 




Figure 4: A+ (solid line) and A~ (point dashed line) particle distributions 
predicted by the R2C model in p—p interactions. Dashed line is the contribu- 
tion to A+ and A~ production coming from the independent fragmentation. 
Experimental data are from Ref. [3]. 

is expected to be small, can only give a marginal contribution, and only at 
high x F (x F — > 1) values, to the Af production in proton-proton interactions. 
Thus, the recombination mechanism seems to be the principal contribution 
to the hadronization process, even more important than independent frag- 
mentation of charm quarks. 
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